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abstract  (Cor\ilnu9  on  te\  fi  5«  sido  If  neceaimry  ond  identify  by  block  nambmr) 

atpr  migration  across  the  thin  epoxy  resin  layer  in  model  adhesive  joints  cause^ 
strongly  inhomogeneous  swelling.  The  swelling  Is  the  principal  cause  of  stress 
systems  which  are  being  examined  by  analysis  of  the  pattern  of  Newton's  rings 
formed  between  a  glass  cover  slip  and  an  optical  flat.  Enhancement  of  water 
migration  by  normal  stresses  of  the  order  of  a  kilobar  is  reported.  The 
development  of  a  self-stressed  fracture  mechanics  test  piece,  de.signed  to  inves-j 
tigate  water  uptake  by  adhesives,  is  reported.  Photoelastic  imaging  of  ultrasonjic 
waves  can  be  utilized  to  investigate  defects  in  opaque  solids.  Instead  of 
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usi-ng  a  p/ezce\ectric  detector,  the  reflected  sound  from  a  defect  is  actually 
observf^  in,a  visualizing  block  of  quartz,  after  propagating  from  the  specimen 

to  the  block  across  an  acoustic  coupling.  Experiments  and  a  computer  model  are 
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1,  Introduction 

During  the  first  of  the  three  years  of  tenure  of  this  contract, 
water  naigration  across  the  thin  epoxy  resin  layer  in  model  adhesive 
joints  was  studied  in  detail.  Migra^on  rates  were  measured  by 
monitoring  the  ingress  of  the  boundary  between  unswollen  (dry) 
and  swollen  (wet)  resin  and  were  faster  than  expected  for  Fickian 
diffusion.  With  no  evidence  to  support  the  presence  of  channels 
and  hence  a  capillary  explanation  for  enhanced  migration,  stress 
systems  were  considered  to  see  if  a  contribution  could  be  attributed 
to  them.  Curing  is  achieved  above  Tg,  the  glass  transition  tempera¬ 
ture  for  the  resin.  Shrinkage  that  accompanies  curing  must  there¬ 
fore  be  largely  relieved  by  viscoelastic  flow.  Differential  contraction 
between  metal  and  resin  during  cooling  from  the  cure  temperature 
is  presumably  at  least  partially  accommodated  by  self-adjustment 
of  the  glue  line  thickness,  although  lateral  tensions  no  doubt  remain. 
A  more  important  source  of  stress  is  the  resin  swelling  that  accom¬ 
panies  water  uptake.  This  swelling  was  shown  to  be  strongly  inho¬ 
mogeneous,  and  is  believed  to  be  the  principal  cause  of  stress  fields 
and  has  been  under  close  study  from  the  inception  of  this  project. 

The  development  of  a  self-stressed  fracture  mechanics  test- 
piece  suitable  for  water  uptake  by  adhesives  was  started  during  the 
first  year  and  has  continued  ever  since.  The  original  test-piece 
design  was  discarded  after  the  first  year  and  has  been  replaced 
by  a  more  realistic  geometry.  This  invention  is  the  subject  of  a 
US  patent  application  filed  by  the  US  Air  Force. 

The  need  to  generate  direct  evidence  for  tlie  presence  and 
locations  of  defects  in  our  adhesive  joints,  both  those  manufactured 
for  swelling  studies  and  those  manufactured  for  fracture  mechanics 
tests,  led  us  into  the  field  of  non-destructive  evaluation.  In 


particular  work  was  started  on  the  photoelastic  imaging  of  interference 
attributable  to  defects  that  are  generated  in  ultrasonic  wave  trains  . 
This  followed  naturally  from  previous  photoelastic  imaging  experi¬ 
ments  relating  to  a  project  on  the  physical  mechanisms  responsible 
for  the  weathering  of  epoxy  resins  t^nd  GFR  epoxy  resins  supported 
by  the  US  Army  grant,  DA.  ERO-76-G-068  .  Attention  was  focussed 
on  a  computer  model  which  investigates  the  waves  originating  from 
a  transducer  probe  and  those  waves  which  emanate  from  a  defect  or 
fault.  Experimentation  has  proceeded  in  parallel  with  computation 
although  difficulties  have  been  experienced.  In  particular,  achieving 
results  at  frequencies  in  excess  of  lOMH^  has  proved  extremely 
difficult  because  of  inherent  increased  attenuation,  and  because  resin 
based  specimens  (adhesive  joints  and  film  composites)  tend  to  dissi¬ 
pate  the  probing  stress  waves  rapidly,  as  a  result  of  the  close  proxi¬ 
mity  of  interfaces  which  the  ultrasound  must  needs  negotiate. 

2.  NEWTON'S  RINGS  EXPERIMENT 
2.1.  Experimentation 

In  geometries  such  as  those  provided  by  adhesive  joints  where, 
over  much  of  the  resin  surface,  direct  access  to  the  environment 
is  restricted,  the  concentration  of  diffused  water  and  the  swelling 
associated  with  it  is  strongly  inhomogeneous.  This  inhomogeneity 
is  conveniently  demonstrated  for  adhesive  films  by  using  a  thin 
and  therefore  flexible  circular  glass  cover  slip  bonded  to  a  massive 
(rigid)  adherend  such  as  a  microscope  slide  or  a  slab  of  metal. 

The  experimental  arrangement  is  shown  in  Figure  1.  Light 
from  a  mercury  vapour  discharge  lamp  is  directed  through  an 
interference  filter,  a  collimating  lens,  a  half-silvered  mirror, 
an  optical  flat,  and  the  cover  slip.  Interference  between  incident 
and  reflected  beams  occurs  within  the  variable  thickness  gap  located 
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between  the  cover  slip  and  optical  flat.  The  interference  pattern  is 
photographed  using  the  light  reflected  into  a  35  mm  camera  by  the 
half  silvered  mirror. 

The  pattern  of  interference  fringes  is  determined  by  the 
geometry  of  the  gap  between  cover  slip  and  optical  flat  and,  to  ensure 
that  changes  in  the  pattern  arise  only  from  distortion  of  the  cover 
slip  caused  by  resin  swelling,  it  is  essential  that  each  run  be  carried 
out  without  disturbing  the  specimen/optical  flat  assembly.  All  of 
the  components  are  set  up  on  an  optical  bench.  Figures  2  and  3 
show  general  views  of  the  ambient  temperature  and  high  temperature 
(  <100°C)  apparatus  respectively. 

It  was  stated  in  the  Interim  Scientific  Report,  29th  September 
1978,  that  enhancement  of  fringe  contrast  as  a  consequence  of  mul¬ 
tiple  beam  interference  can  be  obtained  by  fully  silvering  the  upper 
surface  of  the  cover  slip  and  partially  silvering  of  the  lower  surface 
of  the  optical  flat.  A  range  of  reflective  coatings  has  now  been 
evaluated  from  which  it  is  evident  that  a  chromium  layer  is  most 
effective,  having  sufficient  water  resistance  to  withstand  a  full 
length  boiling  water  test.  The  main  disadvantages  of  a  chromium 
layer  are  that  it  obscures  direct  observations  of  the  adhesive  and  is 
difficult  to  remove  should  it  need  renewing.  An  alternative  and  very 
Successful  method  for  securing  multiple  beam  interference  and  hence 
for  improving  fringe  definition  has  been  realised  by  using  a  wedge 
shaped  optical  flat.  With  a  wedge  angle  of  4°  the  image  of  the 
Newton's  rings  does  not  coincide  with  reflections  from  other  glass 
or  water  surfaces.  The  wedge  is  made  from  flint  glass  with  refrac¬ 
tive  index  1.  65  and  has  been  incorporated  into  a  mark  II  version  of 
the  hot  water  apparatus  designed  to  permit  undisturbed  tests  at  any 
fixed  temperature  between  ambient  and  100°C.  This  is  shown  in 


4. 


Figure  4. 

2.  2  Results  and  Discussion. 

The  first  long  term  test  was  made  at  room  temperature  on  a 
joint  between  a  cover  slip  and  cleaned,  but  not  anodised.  SIC 
aluminium,  A  general  purpose  epoxy  resin  containing  entrapped  air 
bubbles  was  used  as  adhesive.  Hereafter,  this  adhesive  will  be 
referred  to  as  adhesive  A.  Its  full  identity,  including  hardener  and 
accelerator,  and  the  curing  schedule  used  for  it  are  given  in  Table  1. 

Later  tests  w’ere  carried  out  on  joints  between  a  cover  slip 
and  substrates  of  both  stainless  steel  and  anodized  titanium  using 
the  commercially  available  adhesive  film  system  designated  B  in 
Table  1.  All  metal  surfaces  were  cleaned  and,  where  applicable, 
treated  in  accordance  with  commercially  accepted  practice.  The 
metal  based  specimens  w'^ere  supplied  by  Westland  Helicopters  Ltd 
in  the  form  of  cleaned  and  pre-primed  coupon  specimens. 

For  tests  employing  the  adhesive  film  systems  designated 
C  and  D  in  Table  1  both  adherends  were  glass.  Preparation  of  the 
glass  cover  slip  surfaces  is  less  well  standardised.  The  fundamen¬ 
tal  problem  is  removal  of  oxides  of  low  surface  energy.  Chemical 
attack  by  the  common  oxidising  agents,  e.g.  nitric  acid  and  hydro¬ 
gen  peroxide  certainly  remove  a  large  fraction  of  these  oxides.  So 
too  doe.s  exposure  to  the  oxidising  region  of  a  naked  flame.  Sand 
blasting  in  order  to  roughen  the  surface  and  promote  mechanical 
linkages,  and  application  of  coupling  agents  (silanes)  have  also 
been  recommended.  However,  no  combination  of  these  treatments 
has  yet  produced  a  joint  which,  when  it  eventually  debonds  after 
prolonged  exposure  to  water  in  the  absence  of  externally  applied 
stress,  leaves  fracture  surfaces  which  do  not  include  large  areas 
of  glass/adhesive  interface. 


Table  1 
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In  joints  manufactui'ed  using  liquid  resin,  a  reproducible  glue 
line  thickness  was  achieved  by  simply  allowing  the  cover  slip  to  sink 
under  its  own  weight.  In  joints  made  using  solid  adhesive  film, 
the  glue  line  thickness  was  controlled  by  clamping  during  the  cure  pro¬ 
cess. 

Adjacent  rings  of  the  same  colour  (black  or  white)  are  loci  of 
points  for  w'hich  the  optical  path  length,  in  the  space  between  cover 
slip  and  optical  flat,  differs  by  one  wavelength.  By  the  same  token, 
a  displacement  in  the  pattern  of  Newton's  rings  by  an  amount  equal 
to  one  ring  width  corresponds  to  a  change  in  path  length  equal  to  one 
w'avelength.  By  observing  changes  in  the  number  of  rings  between 
fixed  markers,  such  as  entrapped  air  bubbles,  displacements  normal 
to  the  joint  during  water  uptake  can  be  measured  to  an  accuracy 
of  ^  /4.  If  required,  displacements  which  are  at  least  as  small  as 
^  /lO  can  be  measured  by  superimposition  of  images  in  order  to 
create  Moire  patterns. 

The  water  concentration,  and  hence  the  swelling  associated 
with  water  uptake,  eventually  saturates  and  thereafter  a  shoulder 
separating  fully  saturated  from  less  than  fully  saturated  resin, 
progressively  moves  inwards  from  the  rim  of  the  joint.  This  is 
shown  schematically  in  Figui  e  5c.  Figure  6  shows  a  sequence  of 
images  obtained,  showing  the  development  of  swelling  in  a  joint 
manufactured  between  anodized  titanium  and  a  glass  cover  slip  using 
adhesive  B  and  immersed  in  distilled  water  at  approximately  80°C. 

To  maintain  contact  with  the  outer  ring  of  uniformly  swollen 
adhesive  (Figure  5c),  the  adherends  would  each  need  to  bend  with 
curvature  that  is  opposite  to  that  inside  the  shoulder.  Failure  to 
adopt  such  "S"  wise  bending  manifests  itself  as  an  interfacial  crack, 
the  occurrence  of  which  is  revealed  as  patches  of  bright  contrast  at 
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the  rim  iii  Figure  6. 

The  positions  of  well  defined  points  in  the  patterns  of  Newton's 

rings  may  be  measured  as  functions  of  time.  Figure  7  shows  a  plot 

of  the  shoulder  defining  the  extent  of  water  saturation,  and  of  the 

1 

crack  edge  as  functions  of  for  an  aluminium  joint  manufac¬ 

tured  with  adhesive  A  and  immersed  in  distilled  water  at  20°C. 

Figure  8  shows  a  plot  of  the  shoulder  defining  the  extent  of 
water  ingress  from  tiie  experiment  reported  in  Figure  6. 

Figure  9  shows  a  similar  plot  for  a  specimen  manufactured 

using  adhesive  B  but  v/ith  a  stainless  steel  substrate  and  for  an 

o 

immersion  water  temperature  of  80  C.  From  the  marked  similarity 
in  the  results  for  anodized  titanium  and  stainless  steel  it  is  concluded 
that  tlie  deformation  is  wholly  within  the  adhesive  and  is  independent 
of  the  adherend.  Moire  patterns,  generated  by  superimposition  of 
interference  patterns  before  and  after  swelling  were  used  to  obtain 
the  data  in  Figure  10.  These  specimens  were  manufactured  from 
adhesives  C  and  D  and  measurements  were  made  at  a  range  of 
temperatures  between  ambient  and  lOO^C. 

If  migration  of  water  occurs  at  low  water  concentrations  and  is 
governed  by  Ficks'  laws,  plots  of  migration  distance  as  a  function 

i 

of  (time)^  should  approximate  to  straight  lines.  The  data  points  for 
specimens  manufactured  with  adhesives  A  and  B  certainly  do  not 
approximate  to  straight  lines.  However  those  manufactured  with 
adhesives  C  and  D  and  which  were  studied  using  the  more  precise 
Moire  measuring  technique  do  sviggest  PTckian  behaviour  although 
the  slopes  of  the  graphs  indicate  actixation  energies  that  are  much 
higher  than  normally  expected  for  epoxy  resins. 


3.  SELF  STRESSED  FRACTURE  MECHANICS  TEST-PIECE 
3.  1.  The  plane  elastic  problem. 


The  disc  consists  of  a  circle  of  aluminium  which  in  plane  polar 
coordinates  is  defined  by  the  relations  0  r  ^  a  and 

It  is  stressed  by  removing  a  sector  0^  a,  2^  4  6.^  2  R  and 

making  an  adhesive  joint  when  the  gap  has  been  closed.  Assuming 
that  plane  stress  conditions  are  met,  it  is  necessary  to  calculate  the 
stress  intensity  factor,  K^.  and  the  formation  energy,  W,  for  the 
crack  c  ^  a.  6=0.  This  has  been  done  by  superimposing  the 
solution.s  of  the  resulting  stress  with  those  of  and  W  for  a  traction 
free  circumference  and  a  symmetrical  pressure  loading. 

The  Airy  function  given  by  Michell  (1899),  Eshelby  (1966), 


■)^(r,0} 


2  2 

a  +b  Inr  +  c  r  +  d  r  Inr 
o  o  o  o 


may  be  used  to  find  the  stress  and  displacement  fields: 

a  (r.  6)  =  b  r'^  +  2  c  +  d  (1  +  2  In  r) 
rr  '  o  o  o 

0^0  (r.e)  =  0 

-2 

^oQ(r,6)  =  -br  +2c  +d  {3  +  21nr) 

Do  on 


3.  1 

3.2 
3.  3 

3.  4 


u  {r,d)=E'^|-b  (l+'')r'^+2c  (l-»^)r  +  2d  In  r  -  d  (l+V^rl 

r  ^  o  o  o  o  >>3^  5 


and 

Uy  (r.e)  =  4  E"^  d^rO 

Since  the  solution  is  required  to  satisfy  the  conditions 


3.6 


Mitchell,  J.  H.  1899  Proc.  Lond.  Math.  Soc.  31,  100. 
Eshelby,  J*  D.  1966  Brit.  J.  Appl.  Phys.  17,  1131  “1135. 
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{i)a  (a,8)  =  a  o(a,9)  =  0.  04  S  ^  (ii)  uJr,  2T0  =r  **  and 
rr  ro  ^ 


(iii)  u  and  are  bounded  at  r  =  0,  we  see  that  b  -  O.d  =£4^8?^ 

'  '  I'  O  O  (J 

and  c  =-  (E>y8W)  In  (a  ^/o  ).  Therefore, 


^(r.0) 


E*‘ 


r  In  (r  /  a  /e  ) 


and 


ffQ^r.O)  =_^!L  I  1  +  In  (r  /  a)| 
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It  follows  that  we  may  write  a  (r)  =  p^f  (r/a)  where 

p  =E>y4^  and  f  (p  )  =  1  +  In 

To  find  and  W  tlie  dimensionless  quantities  ^  =  r/a  {already 
introduced)  may  be  used 


o 


u(p.0)=  E  U^(r,' 


Po 


and  V  (^,0)  =  E  Ug(r,8) 


p  (1+V)a 


P  (l+^)a 


o  o 

(hen,  making  use  of  a  symmetry  argument,  it  is  easily  shown  that  the 
problem  to  find  and  W  may  be  stated  as  follows: 

Solve  the  dimensionless,  plane  stress  equations  of  elasticity  in  the 
semi-disc  O^p^l,  subject  to  the  conditions 


(i)  The  stresses  and  displacements  are  bounded  at  r  =  0. 

(ii)  Spp  (1,0)  =  (1,0)  -  0  .  0  ^  0  ^  JT 

m  ®p0(P.O)  =  Spg(p,7r)  =  0,  04  p  ^  1 

(iv)  V  (p,  IT)  =  0  ,  0  ^  p  ^  1 

(v)  V  (p,  0)  =  0  ,  0  ^  p  ^  c 

(Vi)  s^g(P,0)  =  -^  +  log  p).  c  ^  P  ^  1 

(vii)  limit  a  v(p.O)  « 

p  1  -  S  p 


By  slightly  modifying  the  results  in  section  2  of  Tw'eed  and  Rooke  (1973) 
it  can  be  shown  that  there  is  a  solution  of  the  equations  of  elasticity,  in  the 
given  r^ion,  which  satisfies  the  conditions  (i)  through  (v)  and  is  such  that 


P(t)  M(P.t)dt 
/[(l-t)(t-c)] 


0  ^  p  4  1 


3.10 


and 
V  (P.O) 


-2_  t  _ P(t)cit 

1  +  J  /[(l-t)(t-c)J 

c 


.  C  ^  P  ^  I 


3.11 


where  P(t)  is  an  arbitrary  function  of  t 
and 


M  (p,t)  = 


_1 _ 

t  -  P 


2  2 
+  (t  -1) 


t(l  -  Pt)‘ 


(1  -  Pt)^ 


(1  -  Pt) 
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1 1  follows  that  conditions  (vi)  and  (vii)  will  be  satisfied  also  if  P(t)  is  a 
solution  of  the  singular  integral  equation 


Tweed  J.  and  Rooke  D.P.  1973  Int.  J.  Engng.  Sci.  11,  65-73. 
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1  /  P(t)M(at)dt 

2  /  ~ — —  =-(l  +  logp),  c  <  p  <  1  3.13 

'  J  /[(i-wt-')] 

c 

with  subsidiary  condition 

P(l)-  0  3.14 


The  stress  intensity  factor  and  the  crack  formation  energy  W 

are  defined  by  the  equations 

, -  F,  (’^.0) 

K  =  limit  J  2  tt  (r  -  ac)  -  - 

I  .  ^  2  A  r 


r->  ac+ 

and  aA 

W  = 


3.15 


/  ["l+logir/a^  u^  (r,  0)  dr 

‘'ac 


3.16 


respectively.  Therefore,  by  3.10  and  3.11,  we  see  that 


K. 


and 


W 


A  2  2 

2  p  a 
o 


t  P(t)  log  t 

y[(l-t)(t-c)j 


dt 


3.17 


3.18 


Let  and  be  the  stress  intensity  factor  and  formation  energy 

respectively  of  a  crack  of  length  2  (a  -  ac)  in  an  infinite  elastic  solid,  when 

the  crack  is  opened  by  a  uniform  pressure  p  ,  then  after  Sneddon  and 

0 

Lowengrub  (1969) 


Sneddon  I.  N.  and  Lowengrub  M.  1969  Crack  Problems  in  the  Classical 
Theory  of  Elasticity  Wiley  J.  and  Sons  N.  Y. 


and 
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Therefore 

K. 
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and 
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’'W  dt 


7f(l-C)  j  /pl-t)  (t-C)J 


Also  of  interest  is  the  strain  energy  release  rate 


(£)  = 


a  w 


E 


where  ^  =  a(l-c)  is  the  crack  length. 
If  we  define 


then 


K. 
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(1-C) 


3.20 

3.21 

3.22 

(jg)  which  is  defined  by 

3.23 

3.24 


o 


(c) 


3.25 


Numerical  solutions  to  equation  3.5  are  found  by  the  Gauss-Chehyshev 

quadrature  technique,  Erdogan  and  Gupta  (1972),  and  the  ratio  KJK  using 

I  o 

Chebyshev-Lagrange  interpolation.  The  Gaussian  quadrature  formula  is 

used  to  determine  W/W  . 

o 

3. 2  Experimentation 

Preliminary  experiments  were  carried  out  using  3  '{7, 6  cm)  diameter 
SIC  aluminium  discs,  1/16"  (0.16  cm)  thick  with  a  5°  sector  cut  out  and  a 
I"  (1.2’  cm)  diameter  central  hole.  The  hole  is  helpful  during  manufacture 
and  reduces  pinching  on  disc  closure.  However,  it  was  found  teat,  taking 
typical  values  for  the  physical  properties  of  NSI  (C)  aluminium,  significant 
areas  of  the  disc  are  stressed  to  levels  above  the  elastic  limit  and  even 
above  the  UTS.  This  serious  objection  to  the  initial  design  was  aggravated 
by  a  buckling  factor.  Mitchell  and  Head  (1961)  calculated  that  a  disclinated 
plate  will  buckle  if  the  angle  exceeds,  i  2buckle  ^  13  (thickness/radius) . 
For  the  original  disc,  /I  =  1°20*  and  would  deform  to  a  cone. 

The  thickness  of  the  disc  was  increased  to  3/16"  (0.48  cm)  which 
ensures  that  closure  after  removal  of  a  5®  sector  can  be  achieved  without 
buckling.  This  decision  to  increase  the  thickness  rather  than  reduce  the 
sector  angle  was  taken  because  of  difficulties  associated  with  manufacturing 
and  gluing  the  disc.  With  the  5^^  sector  cut  out,  gap  closure  with  a  custom 
built  chain  wrench,  was  found  to  cause  the  centre  to  close  too  quickly. 
Parallel  sided  gaps  were  examined  with  gaps  varying  from  0.098"  (2.5  mm) 
to  0,019"  (0.5  mm).  With  gap  sizes  greater  than  0.046"  (1.15  mm)  the 
rims  of  the  disc  pinched  thereby  preventing  closure.  Below  0.  039"  (1  mm) 
the  pinching  effect  while  still  present,  was  less  noticeable.  The  most 

Erdogan  F.  and  Gupta,  G.  D.  1972  Q.  Appl.  Math.  29,  523. 

Mitchell  and  Head,  (1961)  J.  Mech.  Phys.  Solids.  9,  131. 


uniform  closure  was  achieved  with  a  wedge  shaped  cut  out-  The  dimensions 
were  calculated  from  empirical  data  collected  from  repeated  opening  and 
closing  of  the  0.039"  <1  mm)  and  0.019"  (.5  mm)  parallel  sided  gaps.  The 
optimum  wedge  dimensions  were  found  to  be  0. 046"  {1 . 2  mm)  at  the 
specimen  periphery  reducing  to  0.032"  (0.8  mm)  at  the  central  hole.  Some 
plastic  deformation  was  observed  to  take  place  on  closure,  thus  demonstrating 
that  the  maximum  amount  of  elastic  energy  was  being  stored  in  the  discs, 
as  required. 

A  thorough  cleansing  schedule  was  devised  to  remove  all  atmospheric 
corrosion  and  included  anodizing  in  phosphoric  acid.  Adhesive  joints  were 
made  within  two  hours  and  therefore  before  the  onset  of  further  attack  by 
the  atmosphere. 

To  date  three  proprietary  adhesive  resin  films  have  been  used 
FM73M,  BSL  312/5,  and  DLS  468/DLS  476.  Two  restraining  plates  were 
used  to  maintain  discs  flat  during  the  closing  operation  and,  to  obtain 
standardisation,  a  torque  wrench  set  at  16  ft.  -lb.  was  employed.  On 
closure,  the  restraining  plates  were  removed  and  the  clamped  disc  was 
placed  in  an  oven  for  curing.  Strict  adherence  to  the  manufacturer's  data 
sheets  was  maintained  for  the  cure  times  and  temperatures.  The  disc  was 
allowed  to  cool  in  a  desiccator  before  releasing  the  clamp  with  the  disc 
once  again  sandwiched  between  the  restraining  plates. 

Starter  cracks  were  not  artificially  introduced.  Instead,  each 
specimen  was  boiled  in  distilled  water  until,  as  a  consequence  of  resin 
degradation,  a  crack  had  nucl^ted.  Thereafter,  the  crack  length  was 
measured  as  a  function  of  time  of  immersion  in  boiling  water.  Since  the 
crack  length  measured  each  side  of  the  disc  invariably  gave  two  values,  the 


average  measurement  was  taken  as  the  true  value.  Attempts  to  introduce 
starter  cracks  by  making  a  small  cut  or  by  immersing  in  liquid  nitrogen 
in  order  to  shrink  the  aluminium  and  hence  increase  the  stress  field,  proved 
unsatisfactory. 

A  few  experiments  were  performed  using  a  perspex  disc  of  similar 
dimensions  to  the  aluminium  specimens.  The  object  here  was  to  use  a 
photoeiastic  method  to  study  the  stress  field  during  crack  growth.  Crack 
initiation  proved  difficult  and  satisfactory  monitoring  of  the  stress  field 
during  propagation  was  not  achieved. 


3.3  Results  and  Discussion 

On  release  of  the  chain-wrench  clamp  from  around  the  stressed  disc, 
crack  lengths  vary  according  to  the  epojgr  resin  adhesive  formulation  used 
to  make  the  joint.  The  crack  always  propagates  the  full  glue-line  length 
(jt/a  =  1)  for  resin  Redux  312/5,  whereas  for  FM  73M  no  crack  is  present 
in  the  joint  (jt/a  =  0).  The  i/a  values  obtained  for  DLS  4G8/DLS  476 
are  0,  0,  0,  0.40  and  0.82. 

Correlation  of  crack  length  with  independently  measured  values  of 
is  required  in  order  to  generate  calibration  curves  for  a  given 


geometry  of  self  stressed  disc.  Available  data  for  the  dry  adhesives 

^  2  2 

show  variations,  e.g.  FM73M  has  values  of  2.  38  kJ/M  and3.67kJ/M  for 
two  different  batches,  Redux  312/5  and  DI.S  468/DLS  476  both  have  ^ 


values  of  2.8  kJ/M  .  Using  as  calibration  points,  the 


I 


IC 


IC 

values  for  dry 


FM73M  and  Redxix  312/5  and  the  stable  crack  lengths  measured  after 

mechanically  introducing  starter  cracks,  two  values  of  i/a  measured  for 

DLS  468/DLS  476  indicate  that  (  fi.  i  is  3.65  +  0, 55  kJ/M^. 

3  IC  max  - 

Conflict  evidently  exists  between  the  published  data  and  the  results 


obtained  here,  with  DLS  4f^/DI,S  476  our  self-stressed  discs  crack 
partially  or  not  at  all.  With  Redux  312/5,  the  bonds  crack  over  the  full 


length  of  the  glue  line.  With  FM73M  no  cracks  are  found  suggesting  that 


IC 


is  greater  than  2.8  kJ/M  .  With  a  view  to  clearing  up  this  dilemma. 


independent  measurements 


ot§ 


IC 


on  the  particular  batches  of  adhesive 


films  used  here  are  being  made  at  the  University  of  Cambridge, 

Results  presenting  the  variation  of  expressed  as  a  fraction  of 


( 


)  ,  with  water  uptake  by  FM73M  are  shown  in  Figures  11  to  14. 


The  specimens  were  immersed  in  boiling  distilled  water.  The  time  for  the 
starter  crack  to  appear  varied  from  5  to  130  mins.  Examination  of  the 
joints  after  total  failure  show  that  propagation  of  the  cracks  occur  through 
the  resin  and  definitely  not  at  the  resin/aluminium  interface. 


4.  PHOTOELASTIC  IMAGING  OF  DEFECTS 

Andrews  and  Wallis  (1977)  were  first  to  demonstrate  that  a  light 
emitting  diode  (LED)  can  be  used  to  illuminate  the  stress  waves  travelling 
in  a  visualizing  block  when  viewed  through  crossed  polars.  The  method  was 
applied  to  the  detection  of  defects  in  opaque  solids  by  Hardy,  Turner  and 
Ashbee  (1978).  There  is  a  strong  reflected  component,  mode  converted, 
from  a  defect  which  can  be  observed  some  distance  from  the  flaw,  after  the 
stress  wave  has  travelled  from  the  specimen  into  the  transparent  block. 

This  v/ave  interacts  with  direct  waves  from  the  transducer  to  give 
interference  patterns  for  which  a  computer  model  has  been  devised  by 
Kitson,  Low  and  Turner  (1979).  All  the  longitudinal  and  shear  waves 
emanating  from  the  transducer,  the  interface  and  the  defect,  with  and 
without  mode  conversion,  are  investigated. 


Andrews  and  Wallis,  1977,  J.  Phys.  E.  10,  95. 

Hardy,  Turner  and  Ashbee,  1978,  Metal  Science,  12,  406. 

Kitson,  Low  and  Turner  (1979)  'Ultrasonics  International  79'  Graz,  Austria. 


4. 1  Computer  Model 

Figure  15  shows  the  representation  of  the  model,  with  a  series  of 
radiating  point  sources  depicting  the  transducer.  An  array  of  points  P(x,y) 
forms  the  visualized  zone  and  was  originally  set  at  31  x  31  to  conform  with 
computer  graphics.  It  is  assumed  that  the  defect  is  perfectly  round,  that 
longitudinal  (L)  and  shear  (S)  waves  emanate  from  the  transducer  point 
sources,  and  that  on  reflection  and  refraction  both  mode  converted  and 
unconverted  waves  emerge.  The  following  relationships  hold: 


sin  g  _  sin  a  sin  y  sin  6 


where  V  and  V  are  the  velocities  in  the  specimen  and  U_  the  shear 

Xj  o  o 

velocity  in  the  visualizing  quartz.  Note  that  the  subscripts  are  interchange¬ 
able  for  the  incident  shear  wave  reflected  as  longitudinal  and  for  a  refracted 
longitudinal  wave.  Relative  intensities  may  be  evaluated  using  the  Zoeppritz 
equations  and  an  attenuation  factor  of  \// x  is  included  for  the  wave 
amplitude  as  it  travels  away  from  the  source.  Hence  we  may  find  the  exact 
path  which  an>  particular  ray  travels  together  with  its  amplitude  at  any  point. 
Each  flash  of  the  LED  is  represented  by  a  sine  wave,  singly  for  single  flash 
and  several  cycles  for  multi-flash. 

The  complete  picture  is  built  up  as  follows:  Each  type  of  wave 
(longitudinal/shcar  (LS),  long/long/shear  (LLS),  etc.)  is  considered  in  turn 
for  each  transducer  point  source.  The  program  calculates  (i)  the  path 
between  the  transducer  point  and  a  point  in  the  rectangular  visualized  array 
(subject  to  the  path  being  possible,  (ii)  when  the  sound  pulse  arrives  at  the 
visualizing  point,  and  (iii)  the  amplitude.  The  amplitude  is  then  multiplied 


by  the  value  of  the  illumination  W'avefoi'm  at  that  specific  time,  and  repeated 
for  all  the  points  in  the  rectangle.  Results  for  all  the  transducer  elements 
are  superimposed,  and  finally  the  whole  process  Treated  for  each  wav^ype 
in  order  to  simulate  the  complete  mterfermce  pattern.  The  final  array  is 
normalized  to  give  the  negative  or  compressive  regions  the  lighter  shades  of 
the  computer  simulation  and  the  positive  or  tensile  regions  the  darker  areas; 
a  normalization  analogous  to  the  quarter  wave  plate. 

Two  main  calculation  procedures  are  used  in  the  program,  REFPATH 
and  DIRCALC.  REFPATH,  the  flow  diagram  is  shown  in  Figure  16,  deals 
with  the  patterns  resulting  from  wavetypes  which  are  reflected  from  the 
defect,  whilst  DIRCALC  handles  those  direct  from  the  transducer.  Any  wave- 
types  can  be  included,  but  it  is  found  experimentally  that  many  are  of  little 
importance  and  only  LLS,  LS,  and  SSS  need  to  be  considered.  It  should  be 
remembered  here  that  the  experimental  arrangement  of  the  photoelastic 
method  helps  enormously  to  reduce  the  ntunber  of  wavetypes  which  contribute 
to  the  visualization.  In  the  direction  of  principal  stress,  i.e.  the  direction 
of  the  stress  wave  propagation,  45^  to  the  polar  axes,  only  longitudinal  waves 
are  visualized,  and  at  45°  to  this  direction,  along  the  polar  axes,  only 
transverse  waves  can  be  seen.  It  is  in  this  direction,  along  the  polar  axes, 
that  we  are  primarily  concerned  as  it  is  here  that  most  information  is 
displayed. 

The  first  attempts  to  display  the  visualized  array  by  computer  graphics 
* 

took  advantage  of  SYMAP  and  then  SYMVU  was  employed  to  give  a  topo¬ 
graphical  representation  of  the  data.  The  former  program  had  the  advantage 
that  direct  comparisons  could  be  made  of  the  photographic  results,  but 
suffered  from  having  only  ten  levels  of  shading,  whereas  the  latter  shows  the 

*SYMAP  is  a  computer  program  for  producing  maps  which  graphically  depict 
spacially  disposed  quantitative  and  qualitative  information.  SYMAP  was 
designed  and  developed  by  the  Laboratory  for  Computer  Graphics  at 
Harvard  University,  U.S.  A. 


dislocations  very  well.  Figure  17  sfaows  the  contour  map  where  the  maximum 
pressure  is  represented  by  a  height  of  3",  with  a  minimum  of  0".  The 
viewing  point  is  at  an  altitude  of  30"  with  an  azimuth  of  335°.  Dislocations  in 
this  context  refer  to  the  line  of  discontinuities  in  wavefronts  which  line  runs 
in  the  direction  of  the  stress  propagation-  A  new  program  was  written 
therefore,  following  the  principles  of  SYMAP,  but  with  facility  for  100 
shadings.  Some  time  was  spent  experimenting  with  differait  overprinted 
characters  in  order  to  get  a  set  of  evenly  varying  shades.  It  was  found  that 
a  subroutine  designed  to  interpolate  the  brightness  between  the  array  points 
gave  a  smoother  picture  when  the  visualization  density  map  was  plotted  on 
the  line  printer. 

4.2  Visualized  sound  experiments 

Figure  18  shows  a  schematic  diagram  of  the  apparatus.  The  LED 
sends  a  single  pulse  (single  flash)  or  flashes  a  number  of  times  in  synchro¬ 
nisation  with  the  ulti-asonic  frequency  (multi-flash)  to  illuminate  the  stress 
field  in  the  quartz  visualizing  block.  The  birefringence  exhibited  by  quartz 
when  it  is  stressed,  results  in  a  small  component  of  light  being  propagated 
perpendicular  to  the  axis  of  the  polarizer  to  give  elliptically  polarised  light. 
This  contrasts  with  transmission  of  the  illuminating  pulses  through  unstressed 
r^ions  of  the  quartz  since  here  the  light  remaimplane  polarised.  The 
analyser,  set  at  90°  to  the  polarizer,  permits  transmission  of  the  small 
component  only  so  that  the  stressed  regions  are  revealed  as  bright  fringes. 
Bright  fringes  arise  in  both  tensile  and  compressive  r^ions,  and,  by  suitably 
setting  the  quarter  wave  plate,  all  the  light  from  the  most  highly  compressed 
region  can  be  eliminated.  Thus  some  light  can  be  seen  for  unstressed  regions 


with  the  tensiie  areas  brighter  still,  the  whole  picture  being  built  up  strobo- 
scopically  at  repetition  rates  up  to  3  kHz.  Full  details  of  the  experimental 
approach  are  given  in  Hardy,  Turner  and  Ashbee  (1978). 

The  experiments  which  have  been  carried  out  in  parallel  with  the 
computer  modelling  have  been  mostly  sii^le  flash,  i.e.  a  single  LED  pulse 
for  illumination.  Single  flash  images  are  difficult  to  record  because  they  are 
so  dim,  but  it  is  necessary  to  identify  the  different  wavetypes.  A  series  of 
experiments  were  conducted  with  no  specimen  in  place  in  order  to  investigate 
the  transducer  for  accurate  computer  modelling.  For  these,  three  flashes 
were  used.  As  the  shear  waves  travel  more  slowly  than  the  longitudinal 
waves,  the  shear  fringes  are  closer  together:  a  fact  which  with  their 
position  simplifies  the  identification  of  any  particular  wave.  To  investigate 
the  transducer  beam,  a  series  of  photographs  was  taken  varying  the  polar 
axes,  thus  compensating  for  the  angtilar  sensitivity  of  the  optical  arrangement. 
The  angular  spread  of  the  beam  was  found  to  be  approximately  90°. 

A  number  of  different  transducers  have  been  tried  including  two 
commercial  designs,  but  the  best  results  have  been  obtained  with  a  transducer 
designed  and  built  in  the  laboratory.  It  consists  essentially  of  a  brass  plate 
0.9"  (2.2  cm)  diameter  and  5X/4  thick  on  to  which  a  0.4"  (1.0  cm)  diameter 
piezo-electric  disc  is  soldered.  However,  the  image  is  still  complicated  by 
secondary  waves  arising  from  reflections  within  the  brass  and  disc. 

After  the  initial  experiments  it  appeared  that  two  wavetypes  were 
contributing  to  interference  patterns  seen  in  the  visualized  region.  Firstly, 
the  wave  reflected  from  the  defect.  This  is  mode  converted  but  crosses  the 
interface  between  specimen  and  visualizing  block  without  a  further  conversion 
(  LSS).  Secondly,  the  direct  wave  originating  as  shear  at  the  transducer  and 


again  crossing  the  interface  without  mode  conversion  (SS).  The  esqjerinients 
chow,  however,  that  the  prominent  shear  waves  leave  the  specimen/qiiartz 
interface  at  roughly  the  same  time  as  the  primary  longitudinals.  To  investi¬ 
gate  this  further,  the  positions  of  the  primary  longitudinal  and  shear  waves 
were  plotted  for  a  series  of  different  delay  times.  Different  specimen  thick¬ 
nesses  were  used  for  the  measurements  as  well  as  readings  taken  without  a 
specimen.  The  results  shown  in  Figure  19  strongly  suggest  that  the  shear 
waves  arise  by  way  of  mode  conversion  at  the  interface,  having  crossed  the 
specimen  as  longitudinal  waves.  Shear  waves  generated  from  the  transducer 
are  apparently  \inimportant.  Different  effects  of  the  quarter  wave  plate  on 
the  two  wavetypes  account  for  the  incorrect  positioning  of  the  specimen/quartz 
interface  in  that  Figure. 

To  establish  a  baseline  from  which  to  judge  images  from  resin 

specimens  a  rar^e  of  steel  specimens  have  been  produced  with  thicknesses 
1  ” 

up  to  1  /4  (31.2  mm)  each  containing  a  transverse  hole  'defect*  having 
diameters  in  the  range  0.04"  (1  mm)  to  0.4”  (10  mm).  A  35  mm  Nikon  F2 
photomic  camera  is  used  to  record  the  results  with  an  eight  minute  exposure 
on  ASA  400  film  for  single  flash  at  different  delay  times. 

The  delay  time  is  measured  from  the  original  ultrasonic  pulse  to  the 
LED  flash,  and  is  the  time  taken  for  the  sound  to  travel  from  the  transducer 
to  the  visualized  position.  The  TTL  control  enables  this  to  be  accurately 
measured,  since  the  delay  monostable  and  thyrister  fired  ultrasonic  wave  are 
triggered  by  the  same  pulse,  A  digital  readout  of  the  delay  is  used  which  can 
be  checked  against  the  time  base  of  the  oscilloscope. 


4 .  3  Results  and  Discussion 

The  results  are  presaited  in  a  graphical  form  which  show  both 
experimental  position  and  the  computer  simulation.  Figures  20  and  21  show 
the  results  for  two  delays,  8.8  and  13.4  ixs  respectively  for  steel  specimens. 
The  dotted  lines  give  the  position  of  the  computer  prediction  and  the  hard 
lines  have  been  traced  from  the  photograph.  They  show  that  the  program 
correctly  predicts  the  wavefront  spread  for  both  the  LSS  and  LS  waves.  It  is 
the  former  which  represents  ’information'  from  the  defect.  Differeices  in 
the  lengths  of  the  wavefronts  can  be  attributed  to  the  differences  in  the 
transmission  coefficients.  The  latter  are  complicated  by  the  acoustic 
coupling  at  the  specimen/quartz  interface  by  a  thin  film  of  water-  Broadening 
of  the  lower  LS  wave  is  caused  by  the  superposition  of  the  LLS  wave  i.e.  the 
reflection  from  the  hole  without  mode  conversion. 

The  experimental  delay  times  against  the  theoretical  time  required  to 
align  the  LS  waves  is  shown  in  Figure  22.  The  gradient  of  unity  within  3% 
demonstrates  further  that  the  program  correctly  predicts  the  LS  position. 

A  sequence  of  tests  was  undertaken  with  thicker  specimens, 

1  " 

1  A  (31.2  mm)  containmg  a  hole  0. 3"  (8  mm)  diameter,  where  the  distance 
between  the  hole  centre  and  the  transducer,  A,  was  varied.  Figures  23,  24 
and  25  give  the  results  for  A  =  0.27"  (7  mm),  0.61"  (15.6  mm)  and 
0.95"  (24.2  mm)  respectively. 

The  above  results  refer  to  single  flash  which  was  regarded  as  the 
best  approach  for  wave  identification,  whereas  the  original  approach  b^an 
with  multiflash.  More  information  is  carried  in  the  multiflash  pictures  as  the 
close  up  of  the  Interference  patterns  from  a  0,04"  (1  mm)  radius  'defect' 
shows  in  Figure  26.  The  multiflash  simulation  of  Figure  27  illustrates  the 


interaction  of  the  LS  wave  and  the  LSS  wave  reflected  from  a  hole  of  0.04*' 

(1  mm)  radius. 

In  composite  materials  it  has  beoi  possitie  to  calculate  both 
longitudinal  and  shear  wave  velocities  by  measuring  the  delays  across  a 
known  specimen  thickness.  Velocity  variations  have  been  detected  during 
water  uptake  and  are  attributed  to  resin  plastication  and  to  post  curing. 
Experiments  on  composite  laminates  containing  holes  have  been  less 
successful.  Attenuation  of  the  wavctypes  is  more  marked,  and  the  shear 
wave  rapidly  loses  energy  as  it  proceeds  across  a  series  of  interfaces.  The 
proUem  is  aggravated  by  higher  frequencies,  resolution  of  the  fringes 
being  impossible  at  frequencies  in  excess  of  ~  10  MHz. 


General  view  of  ambient  temperature  apparatus. 
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Figure  4 


Mark  II  hot  water  apparatus. 


Sequence  of  interference  patterns  showing  the  swelling  and  debonding 
film  adhesive  during  immersion  in  water  at  80°C. 


Figure  7  Migration  of  the  shoulder  in  a  Newton's  ring  (□  )  and  of  a  point  on  the  debonding  crack  (O  )• 

20®  water. 


Figure  8  Diffusion  rate  for  an  FM  1000  adhesive  joint  between  anodized  titanium 
and  glass  cover  slip. 


Figure  9,  Diffusion  rate  for  an  FM  1000  adhesive  joint  between 
stainless  steel  and  glass  cover  slip 


Migration  distance 


Figure  12  Plot  of 


Transducer 
1  Centre 


Figure  16 


Flow  diagram  of  program  REFPATII. 


■e  18  Schematic  diagram  of  photoeiastic  system. 


Photograph  - Computer  prediction 

Photograph- Printout  Comparison 


(x,y)  =  (1*50, 2  00)  cm 


Photograph  - Computer  prediction 

Photograph- Printout  Comparison 

Delay  =  13-.4jus 


Figure  21 


Computer  simulation 
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Photograph  - Computer  predict 

Photograph- Printout  Comparison 
A  =  7-0  mm 


Figure  23 


Computer  simulation 


Photograph  - Computer  prediction 

Photograph- Printout  Comparison 

A  =  15-6  mm 


Figure  24 


Computer  simulation 


45<  4  cm 


Photograph  - Computer  prediction 

Photograph- Printout  Comparison 
A  =  24’2  mm 


Figure  25 


Computer  simulation 


